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Abstract. This paper discusses finite-dimensional (Markovian) reaUzations 
(FDRs) for Heath-Jarrow-Morton interest rate models. We consider a d- 
dimensional driving Brownian motion and stochastic volatility structures that 
are non-degenerate smooth functionals of the current forward rate. In a recent 
paper, Bjork and Svensson give sufficient and necessary conditions for the 
existence of FDRs within a particular Hilbert space setup. We extend their 
framework, provide new results on the geometry of the implied FDRs and 
classify all of them. In particular, we prove their conjecture that every short 
rate realization is 2-dimensional. More generally, we show that all generic 
FDRs are at least (d + l)-dimensional and that all generic FDRs are affine. 
As an illustration we sketch an interest rate model, which goes well with the 
Svensson curve-fitting method. These results cannot be obtained in the Bjork— 
Svensson setting. 

A substantial part of this paper is devoted to analysis on Frechet spaces, 
where we derive a Frobenius theorem. Though we only consider stochastic 
equations in the HJM-framework, many of the results carry over to a more 
general setup. 



1. Introduction 

This paper discusses the existence of finite-dimensional forward rate models 
within the framework of Heath, Jarrow and Morton (henceforth HJM) Fol- 
lowing Musiela [|l5j, we denote by rt{x) the instantaneous forward rate at time t 
for maturity t + x. The price at time i of a zero coupon bond with maturity T is 



then given by 



pT-t ^ 
P(i,T) = exp I - / rtix)dx 



It is shown in that essentially any HJM- type model for rt (x) can be considered 
as a stochastic equation 

f drt = [Art + anjAiit, n)) dt + a{t, n) dWt 

on a separable Hilbert space H which is characterized by the following three prop- 
erties: 

(HI) H C C(K>o; M) with continuous embedding; that is, the pointwise evaluations 
evx ■ h t-^ h{x) are continuous linear functionals. 



Date: December 11, 2000 (first draft); February 1, 2008 (this draft). 

^Typically, H consists of equivalence classes of functions, which we shall identify with their 
continuous representatives, respectively. 
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(H2) The family of right-shifts, S{t)f{x) = f{t + x), forms a strongly continuous 

semigroup {S{t) \ t G IR>o} on H. 
(H3) There exists a closed subspace Hq of H such that 

S{f,g){x) -.^fix) r g{v)dri, 
Jo 

defines a continuous bilinear mapping S : Hq x Hq —>■ H. 

We write shortly S{f) for S{f, /). 

As an illustration we shall always have the following example in mind (see 
Section 5]). 

Example 1. Let w : M>o [l,oo) be a non- decreasing C^-function such that 
w^^^^ £ L^(R>o). We may think of w{x) — e°^^ or w{x) — + for a > or 
a > 3, respectively. The space consisting of absolutely continuous functions h 
on M>o and equipped with the norm 

\\h\\l := \h{0)\^+ [ \h'{x)\''w{x)dx 



is a Hilbert space which satisfies (H1)-(H2). Property (H3) is satisfied for Hq — 
Hi {h e H^. I lim.^eo h{x) = 0}. 

Throughout we are given a filtered probability space {Q,T,{Tt),P) satisfying 
the usual conditions. The predictable cr-field is denoted by V. Now we can give 



a meaning to the coefficients in (1.1). Here W = {W^, . . . ^W^), d £ N, denotes 
a standard d-dimensional Brownian motion. The volatility coefficient a{t,uj,h) ~ 
{a^{t, u},h),. . . , (j'^{t, uj, h)) is a measurable mapping from (M>o x f2 x H,'P®B{H)) 
into {H^,B{H^)), well-defining 

d 

aHjM{t,uj,h) := y^^S{a-^t,uj,h)), 

a measurable mapping from (M>o x Vt x H,V ® B{H)) into (iJ, B{H)). 

The operator A is the infinitesimal generator of the semigroup {S{t) \ t G M>o}. 
Its domain in H is denoted by D{A), and we write A* for its adjoint. Both oper- 
ators, A and A*, are closed. It is easy to see that D{A) c {h e H n C^{R>o;M.) \ 
h' e H} and Ah — h' . Without much loss of generality we shall in fact assume 

(H4) D{A) ^ {h e H r\C\R>a;R) \ h' e H}. 

Also (114) is satisfied for the spaces ffu, from Example |l|. 

We distinguish, in decreasing order of generality, between (local) mild, weak 



and strong solutions to equation (1.1) and its time-shifted versions. The reader is 
referred to ||] or Q for the precise definitions. 



Definition 1. A subset V of H is called locally invariant for (1.1) if for all space 



time initial points (io, ^o) G II^>o ^ V ? there exists a continuous local weak solution 



j-{to,ho) iq (1.1) with lifetime r = T{to,ha) such that 



r['°;""^ev, yteR>o. 

For the notion of a finite-dimensional C'^-submanifold of a Hilbert space we refer 
to [pl. We restate §, Theorem 31. 
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Theorem 1. Suppose that a{t,Lu,h) is continuously dijferentiable in h and right- 
continuous in t. Let Ai be an m- dimensional -submanifold of H. Then the 
following conditions are equivalent: 

i) Ai is locally invariant for 

ii) M (ZD{A) and 

1 

^i{t,UJ,h) ■.^Ah + aHjMit,uj,h)--Y,Da\t,uj,h)a^{t,uj,h) E T^M (1.2) 

aHt,uj,h) eThM, j = l,...,d, (1.3) 
for all (i, h) S R>o x M, for F-a.e. lo. 



Hence the stochastic invariance problem to ( |l . l| ) is equivalent to the deterministic 
invariance problems (for any lo apart from a P-nuUset) related to the ^-dependent 
vector fields {/i, cr^, . . . , a'^}. 

Note that Theorem |l| provides conditions for the invariance of a single given 
submanifold A4. It does no t say anything about the existence of finite-dimensional 
invariant submanifolds for (|l.l| ) or {/i, cr^, . . . , a'^}, respectively. This issue will be 
exploited in the present paper. 

We remark that the existence of finite-dimensional invariant submanifolds is 
essentially equivalent to the existence of finite-dimensional realizations (FDR) for 
( Pl| ) in the following sense (see [|[ Theorem 6.4.1]). 

The ore m 2. Let a and Ai be as in Theorem [J. Suppose Ai is locally invariant 



for (LJ). Then, for any space-time initial point {to,ho) G K>o x A4, there exists a 
C^-map (j) : R™ ^ U Ai, where U is an open neighborhood of ho in H, and an 
M."^ -valued diffusion process Z such that 

ril^;""'^ ^ ^{Zt^r), Vt£R>o, (1.4) 
for some stopping time r > 0. 



Following [|3[ we say that equation (1.1) has a local m-dimensional realization at 



{to, ho) if ( |l.4| ) holds. We say that (1.1) generically admits an FDR at {to, ho) if 
( |l.l| ) has a local m-dimensional realization for initial points in a neighborhood of 

{to,ho). 

Remark 1. In view of the preceding remarks we shall suppress the uj-dependence 
from now on, making the convention that all subsequent statements (have to) hold 
simultaneously for all uj (apart from a V-nullset). 

Remark 2. The forward rate is a mathematical idealization rather than an object 
that can directly be observed on the market. For statistical inference one usually 
estimates the current forward curve by fitting a parametrized family Q of curves 
to the finitely given many market data. Such a family Q is generically a finite- 
dimensional immersed submanifold in H . For the sake of consistency of the curve- 



fitting procedure with any given stochastic interest rate model of the form (1.1), it is 
vital that Q is invariant for the flow of ^l.\) . If this is not the case, it is unclear what 
the parameters of the daily fitting should mean in terms of the model. Consequently, 
the choice of the volatility structure, determines the consistent curve- 

fitting procedures. Volatility structures that do not admit invariant .submanifolds 
(i.e. FDRs) should therefore be avoided. 
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The FDR-problem consists of finding sufficient conditions on ct^, . . . , ct'' for the 
existence of generic FDRs. Bjork et al [Q, |^ translated this into an appropriate 
geometric language. In Q they completely solved the FDR-problem for a very 
particular choice of H . Their key argument is the classical Frobenius theorem, 
since they are looking for foliations (which is the appropriate notion for the FDR- 
problem on Hilbert spaces). Therefore they define a Hilbert space, TL, on which 
A — d/dx is a, bounded linear operator. As a consequence Ti consists solely of 
entire analytic functions (see |||, Proposition 4.2]). It is well known however that 
the forward rates implied by a Cox-IngersoU-Ross (CIR) Q short rate model are 
of the form rt{x) = go(a^) + '<'t{^)gi{x) where 

go (a;) = d and Qiix) — -77, 

for some a,6, c > and d > (see e.g. |||, Section 7.4.1]). Since both go and gi, 
when extended to C, have a singularity at a; = (log(c) -I-itt) / a, they cannot be entire 
analytic. Hence the CIR forward rates do not belong to TO. This fact seems to be 
inconsistent with some of the results in |^ (see Remark 7.4 therein). Since CIR- 
type structures appear generically in HJM-framework, the Bjork-Svensson setting 
is too narrow, even though all geometric ideas are already formulated. 

To overcome this difficulty we have to choose a larger forward curve space. How- 
ever, we cannot get around the Frobenius theorem, which requires smoothness of the 
vector fields {/i, cr^, . . . , cr'*}. The problem here is that A is typically an unbounded 
operator (the choice of H in is exactly made to overcome this problem). The 
appropriate framework for an extended version of the Frobenius theorem is given 
by the Frechet space 

D{A°°) := Pi DiA""), 
equipped with the family of seminorms 

n 

P,,{h) = ^"'^Ilff' " e No. 

m=0 

Indeed, D(A°°) is essentially the largest subspace of H which is invariant for A. This 
calls for a Frobenius theorem in Frechet spaces, which is definitively not straight- 
forward, since the inverse function theorem fails in Frechet spaces in general (see 
|T2| for example). Denote by ^0 '■ D{Aq) Hq the restriction of A to Hq- That 
is, D{Aq) = {h e D{A) n Hq I Ah e Hq}. The definition of the Frechet space 
D{A'^) = n„gN^(^o) is obvious. The next result follows immediately from (HI), 
(H3) and (H4). 

Lemma 1. For any f,g E D{Aq) we have S{f,g) G D{A) and 

AS{f, g) = S{Af, g) + S{f, Ag) + / ev^{g). 

Hence S : D(A^) x D{A'^) — > D{A°^) is a continuous bilinear mapping. 

The preceding specifications for a are still too general for concrete implemen- 
tations. We actually have the idea of cr being sensitive with respect to func- 
tionals of the forward curve. That is, a^{t,h) — cj)^ {t,£i{h), . . . ,£q{h)), for some 
q £ N, where 4^ : R>o x — + D{A'^) is a smooth map and ii, . . . ,£q denote 
continuous linear functionals on H (or even on C(R>o;K)). We may think of 
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ii{h) = i^/xi) J^^ h{ri) dr] (benchmark yields) or £i{h) = evxi{h) (benchmark for- 
ward rates). This idea is (generahzed and) expressed in terms of the following 
regularity and non-degeneracy assumptions: 

(Al) a{t,h) = (l>{t,£{h)) with £ : H ^ B smooth linear and </) : B D{A^)'^ 

smooth, so it is a Banach-map (see Definition ^ below). 
(A2) {£, £oA): DiA"^) ^ B x B is a Hnear open map. 

(A3) {0^(t, y), . . . , 'p'^it, y)} are linearly independent, for all {t, y) G K>o x B. 
(A4) A is unbounded, that is, D{A) is a strict subset of H. 

We believe that this setup is flexible enough to capture any reasonable HJM-type 
interest rate model. Combining Lemma ^ and (Al) yields 

Lemma 2. S{a^ {■)) : K>o x H ^ D{A°°) is a Banach-map, for any j ~ 1, . . . ,d, 
hence also ctHjM- 

Assumption (A4) is equivalent to the following property of A. 

Lemma 3. Let A be the generator of a strongly continuous semigroup S on a 
Banach space Y , then the operator A : D{A°°) — s- D{A°^) is a Banach-map if and 
only if A : Y ^ Y is bounded. 

Proof. Let y be a Banach space and A the generator of a strongly continuous 
semigroup S. If ^ : D{A°°) — > D{A°°) is in a neighborhood of C/ a point xq 
a Banach map, then there are smooth mappings R : U C D{A°^) X and 
Q : V C X ^ D{A°°) such that A ^ Q o R and X is a Banach space. By 
differentiation we obtain at xq 

A^ DQ{xo) ■ DR{xa) 

which means in particular by continuity that there exists n > such that DR{xo) 
can be extended continuously to a linear mapping DR{xo) : D{A"-) X. So 
A : D{A") D{A") is a continuous mapping. We recall the Sobolcv Hierarchy 
for strongly continuous semigroups (see p^) defined by the following commutative 
diagram 



Y ^Y 



D{A)^^D{A) 



R{\) 



5f) 



R(\) 



D{A^)—^D{A 



Here R{X) := (A — A)~^ denotes the resolvent at a point of the resolvent set, 
which defines an isomorphism from D{A^) to D(A"+^). The semigroups S*'"^ are 
defined by restriction and are strongly continuous in the respective topologies. The 
generator of S'(") is A restricted to D(^("+i)). If A were continuous on Z?(^("+^^), 
then we could climb up the diagram by the resolvent isomorphisms. A continuous 
on £)(A(")) implies that S'(") is a smooth group, so by climbing up through the 
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isomorphisms S'^^^' is a smooth group and therefore the infinitesimal generator is 
continuous, since it is everywhere defined, by the closed graph theorem. □ 

We will demonstrate that the Frobenius integrability condition implies a very 
particular geometry of the invariant submanifolds - loosely speaking, each of them is 
a band of copies of an affine submanifold. This is similar to Hamiltonian mechanics, 
where the existence of several conservation laws forces the flow to be viable not only 
on a low dimensional submanifold, but even on a special one, namely a torus. 



Remark 3. Although we focus on the particular stochastic equation (1.1), many 
arguments can be carried over to more general stochastic equations in the spirit of 
Da Prato and Zabczyk |^ (see also Remark^^). 

2. Analysis on Frechet Spaces 

A Frechet space i? is a completely metrizable locally convex vector space. On a 
Frechet space E we are equivalently given a sequence (which can be chosen increas- 
ing) of seminorms {pn}n>i generating a complete locally convex topology. For the 
construction of a differential calculus we need the concept of differentiable curves 
into the Frechet space and differentiable mappings on open subsets of a Frechet 
space. To obtain a completely working differential calculus we should enlarge the 
category of vector spaces to so called convenient vector spaces, however, for our 
purposes we can directly give the definitions on Frechet spaces (see for all de- 
tails). We remark that already on Frechet spaces the situation concerning analysis 



was complicated and unclear until convenient calculus was invented (see |12 for 
extensive historical remarks). 

We denote the set of continuous linear functionals hy E'. A subset S of a Frechet 
space is bounded if and only if 1{B) is a bounded subset of M for all I £ E' . A 
multilinear mapping m : i?i x ... x £"„ ^ F is called bounded if bounded sets Bi x 
... X Bn are mapped to bounded subsets of the Frechet space F. On Frechet spaces 
bounded linear mappings are continuous (see or any textbook on locally convex 
spaces for the proof). In general this is not the case, however continuous linear 
functionals are clearly bounded. The set of bounded linear operators on a Frechet 
space E is denoted by L{E). We prefer the notion "bounded" to "continuous" in 
the case of linear operators, since it is better adapted to convenient calculus (the 
following result is true for all convenient spaces). 

Definition 2. Let E be a Frechet space, then c : K — > E' is called smooth if all 
derivatives exist as limits of difference quotients. The set of smooth curves is de- 
noted by C°°(]R, E). A mapping f : U C E R is called smooth if U is open and 
foce C°°(M,R) for all c e C°°(R,E). 



Theorem 3. Let E,G,H be Frechet spaces, U C E, V G G open subsets, then we 
obtain: 

i) A curve c : M. —^ E is smooth if and only if it is weakly smooth, i.e. I o c £ 
C°°(]R,R) for all I £ E' . 

ii) Multilinear mappings are smooth if and only if they are bounded. 

iii) If f : U G is smooth, then df : U x E ^ G is smooth and bounded linear in 
the second component, where 

df{x,v) := ^\t=ofix + tv) 
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iv) The chain rule holds. 

v) Taylor's formula is true, where higher derivatives are defined as usual: 
for all n eN. 

vi) Given c : R — > L{E), then c : R x E ^ E is smooth if and only if eVy o c G 
C°°(M, E) for all v Cz E, where evy : L{E) —^ E denotes the point evaluation. 

vii) A family Oi G L{E) for i G I of bounded linear operators is uniformly bounded, 
i.e. for all bounded subsets B G E there is a bounded subset C d E such that 
for alii I we have ai{B) C C, if and only if it is pointwise bounded ("smooth 
uniform boundedness principle" ) . 

viii) Let B be a Banach space and f : U G be a smooth mapping, then f is locally 
Lipschitz, i.e. for every point x & U there is a neighborhood V <ZU with 

sLLA ^I^l-j; 3^ w G V\ is bounded in G 

\\x-y\\ 

Remark 4. The introduction of analysis on Frechet spaces via definition 2 is not 
limited to Frechet spaces, but works on all locally convex spaces ( so called convenient 
vector spaces) where weakly smooth curves are smooth. The necessity to work on 
these spaces instead of working on Frechet spaces arises from the fact that spaces 
like C'°°{E,F) or even L{E) generically are not Frechet space for Frechet spaces 
E,F. 

Even in the case o/M^ it is not obvious that the smooth functions of definition 2 
are smooth in the classical sense. Notice that convenient calculus completely solves 
the problem, how to do analysis on locally convex vector spaces, and provides a very 
powerful tool in concrete calculations (see for many examples). 

Notice that on Frechet spaces the notion of C"-curves is useful, too, and that we 
are also given the chain rule with smooth functions by the uniform boundedness 
principle: 

mt)) - /(c(o)) _ r ^^(,(0) + _ ,(0))) . 



t Jo t 

^Df{c{0))-c'iO) 

as t ^ 0. Furthermore we obtain the following useful result by Taylor's formula 
and uniform boundedness (see |^ for details). 

Theorem 4. Let E, G be Frechet spaces and f : U C E ^ G a smooth mapping. 
Given xq U and a continuous seminorm q on G, then there is a continuous 
seminorm p on E and £ > such that 

Qifix) - f{y)) < p{x - y) 
if p{x — xq) < e and p{y — xq) < e. 

Concerning differential equations, there are possible counterexamples on non- 
normable Frechet spaces in all directions, which causes some problems in the foun- 
dations of differential geometry (see and the excellent review article Nev- 
ertheless a useful generalization of the existence theorem for differential equations 
on Banach spaces is given by the following Banach mapping principle (see M for 
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details) . Given P : U C E ^ E smooth, we are looking for solutions of the ordinary 
differential equation 

X :] — e,e[^ U smooth 
Ixit) = Pixit)) 

x{0) e u 

If for any initial value y in a small neighborhood V oi xq £ U there is a unique 
smooth solution t i— > Xy{t) for t g] — £,£[ depending smoothly on y, then Fl{t,y) := 
Xy{t) defines a local flow, i.e. a smooth map 

Fl ■]-e,e[xV E 

Fl{0,y) = y 

Fl{t,Fl{s,y))^Fl{s + t,y) 

if s, t,s + t e] ~ £, e[ and Fl{s, y) G V. If there is a local flow around Xq G U (this 
shall mean once and for all: "in an open neighborhood of xq"), then the differential 
equation is uniquely solvable around xq £ U and the dependence on initial values is 
smooth. Let U be connected, then a maximal local flow Fl associated to P is a local 
flow defined on a connected open set W C M. x U such that for any integral curve 
X : I ^ U with a:{0) = y and / connected we have / x {y} G W and Fl{t, y) = x{t) 
for t G I. UP admits a local flow around any point xq G U, then it admits a 
maximal local flow. A maximal local flow satisfies 

{t, Fl{s, y)) eW a and only if (s + t,y) gW 

for all (s, y) gW and i e M. In this case Fl{t, Fl{s, y)) — Fl{s + 1, y) (on flows see 
for example [p2|). 

Definition 3. Given a Frechet space E, a smooth mapping P : U d E ^ E 
is called a Banach map if there are smooth (not necessarily linear!) mappings 
R:U G E ^ X and Q -.V G X E such that P = Q o R. 



U ge 




Theorem 5. We denote by B{U) the set of Banach map vector fields, by X{U) all 
vector fields on an open subset of the Frechet space E. Then B(U) is a C°^'{U,M)- 
submodule ofX(U). 

Proof We have to show that for f,g G C°°{U,R) and Pi,P2 G B{U) the linear 
combination fPi + 5P2 G B{U). Given Pi = Qi o Ri for i = 1,2 with intermediate 
Banach space Xi , then fPi+gP2 = QoR with Q : R'^ XV1XV2 G x XiX X2 ^ E 
and R : U R"^ X Xi X X2 such that 

Q{r, s, vi,V2) = rQi{vi) + 5(52(^2) 
Rix) = if{x),g{x),R,{x),R2ix)) 

So the sum fPi + 5P2 is a Banach map and therefore the set of all Banach map 
vector fields carries the asserted submodule structure. □ 



FINITE DIMENSIONAL REALIZATIONS FOR STOCHASTIC EQUATIONS 



9 



Theorem 6 (Banach mapping principle). Let P : U <Z E E he a Banach map, 
then P admits a maximal local flow. 

Proof. For the proof see □ 

Remark 5. Parameters and time- dependence are treated in the following way. 
Given an open subset of parameters Z C Y of a Banach space and P : I x Z xU 
E, where I is a open set in R and U is open in E, such that Pt^p = Qt,p ° Rt.p, 
where Q and R depend smoothly on time and parameters, P admits a unique smooth 
solution for any initial value Xq E U at any time point tQ (z I depending smoothly 
on parameters, time and initial values. For the proof we look at the extended space 
G -.^RxY xE with P{t,p,x) = (1, 0, P(t,p, x)) and 

Q{t,p,z) = (l,0,QtAx)) 

R{t,p,x) = {t,p,Rt,p{x)) 

with Banach space X :^ R x Y x X . 

We can replace in the above definition the interval ]—£,£[ by [0, e[ to obtain local 
semiflows if and only if the differential equation admits unique solutions around an 
initial value depending smoothly on the initial values. If for all x G U there is a 
local semiflow around x, then there is a maximal local semiflow. A maximal local 
semiflow Fl : W ^ E satisfies 

{t, Fl{s, y)) eW if and only if is + t,y) gW 

for all {s,y) € W and t G R>q. In this case Fl{t, Fl{s,y)) = Fl{s + t,y). The 
notion of a semiflow is redundant on finite-dimensional vector spaces. 

We are in particular interested in special types of differential equations on 
Frechet spaces, namely Banach map perturbated bounded linear equations. Given 
a bounded linear operator A : E ^ E, the abstract Cauchy problem associated to 
A is given by the differential equation associated to A. We assume that there is a 
smooth semigroup of bounded linear operators S : M>o x E ^ E such that 

lim — = Ax 

tio t 

which is a global semiflow for the linear vector field A. Notice that the theory of 
bounded linear operators on Frechet spaces contains as a special case Hille-Yosida- 
Theory of unbounded operators on Banach spaces (see for example Given a 

Banach map P : U C E ^ E with splitting P ~ Q o R we want to investigate the 
solutions of 

^^x{t)=Ax{t)+P{xit)) 

We do this by constructing a solution to the integral equation arising from variation 
of constants 

Fl{t,x) = Stx+ [ St-sP{Fl{s,x))ds 
Jo 

for small positive time intervals and an open neighborhood of a given initial value. 
Given xq & U there exists - due to theorem 4 - a seminorm p on E and 5 > such 
that 

-i?(a;2)|| < p{xi ~ X2) 
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for p{xi — xq) < 6 and i = 1,2, where ||.|| denotes the norm on X. Furthermore 
given yo G X, then for any seminorm q on F there are constants Cq and dq such 
that 

q{Q{yi)-Q{y2))<Cg\\y,-y2\\ 

for Wvi — yo\\ < Sq and i = 1,2. By the uniform boundedness principle the set of 
continuous linear operators {St}o<t<T is uniformly bounded for any fixed T > 0, 
i.e. for any seminorm p on E there is a seminorm qp such that 

p{Stx) < qp{x) 

for t < r. We denote by C([0,e],X) continuous curves on the interval [0,£] to X, 
yo ■— R{xo)- Without any restriction we can assume that xq = and yo = by 
translations. We can then define a mapping 

M:U' cExV C C{[0,e],X) V 

such that M{x,g)(t) = R{Stx + St-sQ{gis))ds) for t € [0,e]. Given g e 
C{[0,e],X) such that \\g{t)\\ < 9 iov < t < e with {g\supt\\g{t)\\ < 9} C V, 
we have 

p{Stx + I St-sQ{9{s))ds) < qp{x) + e{qp{Q{y„)) + Cq^9) 
Jo 

provided 9 < Sq^. This can be made smaller than ^ if e is appropriately small and 
U' := {x £ E with qp{x) < r]} with r] appropriately small. In particular Cq^e < 1. 
If we assume these conditions, then M is well defined, continuous and furthermore 

snp\\M{x,gi){t) - M{x,g2m\\ < £s^Pqp{Q{9i{t)) - Q{92{t))) < 
t t 

<Cq^s sup \\gi{t)-g2{t)\\ 
t 

Consequently M{x, .) is a contraction in V with contraction constant bounded 
uniformly in x € ?7' by a constant strictly smaller than 1. It follows that there is a 
unique g{t,x) for any x €U' depending continuously on x, such that 

M{x,g)=g 

by the contraction mapping theorem. We define 

Fl{t, x) := Stx + I St-sQ{9{s, x))ds 
Jo 

and obtain 

Fl{t,x) = Stx+ [ St-sP{Fl{s,x))ds 
Jo 

since R{Fl{t,x)) = g{t,x) by construction. By induction, smoothness with respect 
to time is easily established and uniqueness follows from the contraction mapping 
theorem. In the time dependent case - with respect to the bounded linear map 
and the Banach map - we can argue in the same manner and obtain a unique 
time-dependent semiflow, smooth with respect to time. 

Concerning smoothness with respect to the initial value, we proceed in the fol- 
lowing way: First we show that there exist directional derivatives: By Taylor's 
formula we obtain 

P{xi) - P{X2) = L{xi,xo) ■ {Xi - X2) 
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where L{xi,X2) ■ h DP{xq + s{xi — Xq)) ■ hds is a Banach map in all three 

variables. So we can solve the system given by 

{xo,xi,h) ^ {Axo + P{xo),Axi + P{xi), Ah + L(xi,X2) ■ h) 

with a flow Fl{t, xq, xi, h) = {Fl{t, xq), Fl{t, xi), M{t, xq, xi, h)), smooth in time 
and continuous in initial values, where the dependence on h is homogenous, so the 
flow can be defined everywhere in h. By uniqueness the identity 

j^iFl{t,xo)~Fl{t,xi)) 

= A{Fl{t, xo) - Fl{t, xi)) + L{Fl{t, xo),Fl{t, xi)) ■ {Fl{t, xq) - Fl{t, xi)) 
leads to 

M{t, Xq, Xi,Xo — Xi) = Fl{t, Xq) — Fl{t, Xi). 

By homogenity in h we obtain the existence of the directional derivatives and its 
continuity in point and direction at the domain of definition. The argument follows 
an argument given in joj. If the dependence were smooth, the derivative with 
respect to the initial value would satisfy the integral equations 

D2Fl{t, x)-v^Stv+ [ St-s ■ DP{Fl{s, x)) ■ D2Fl{t, x) ■ vds 
Jo 

Fl{t,x) ^ Stx+ I St-sP{Fl{s,x))ds 
Jo 

providing an equation of the same type with Banach mapping 

ix,y)^{P{x),DPix)-y) 

and semigroup St © St for t > 0. This equation has by the above procedure a 
unique solution, continuous with respect to initial values. If we integrate u t— » 
D2Fl(t,x + uv) ■ V with respect to u we obtain Fl{t,x + uv) — Flit, x) by uniqueness 
and Taylor's formula for P. Thus we can conclude by Theorem 12.8. of |Q and 
induction. There is one argument hidden, namely, that the domain in t might 
shrink. We overcome this difficulty by the local existence result for the directional 
derivatives, which means that we can construct solutions for the derivative-equation 
by taking original solutions. 

In an analogue manner we proceed in the time dependent case. Consequently we 
obtain smooth semiflows in all possible cases, which proves the following theorem: 

Theorem 7. Let E he a Frechet space and I an interval, A : I x E E a smooth 
curve of bounded linear operators, such that there is a local time dependent semiflow 
Ev(t, s, x) for t > s > in appropriate domains and x G E satisfying 

-^Ev{t, s, x) = A{t)Ev{t, s, x) 
Ev{t, t,x) = x 

for t > s > and x £ E. If P : I x U ^ E is a time- dependent family of Banach 
mappings, then there is a time dependent local semiflow Fl(t,s,x) for t > s > 
appropriate and x £ E satisfying 

—Flit, s, x) = A{t)Fl{t, s, x) + Pit, Flit, s, x)) 
dt 

Fl{t,t,x) = X 
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for t > s > and x E E. 

Remark 6. Parameter dependence can be easily treated by the above methods on 
the extended phase space. 

3. SUBMANIFOLDS AND WEAK FOLIATIONS 

We are interested in the geometry generated by a finite number of vector fields 
given on an open subset of a Frecfiet space E. Therefore we need the notions of 
finite-dimensional submanifolds (with boundary) of a Frechet space (see for all 
details and more). First we introduce the notion of a smooth Frechet manifold. 

A chart on a set M is a bijective mapping u : U ^ 'u-iU) C Ejj, where Ejj 
is a Frechet space and U C M, u{U) C Eu is open. We shall denote a chart 
by (t/, m) or (u,w(t/)). For two charts {Ua,Ua), {Up, up) the chart changing are 
given by u^p := o : upiUap) "^(J/q/j), where := Ua n Up. An 
atlas is a collection of charts such that the Ua form a cover of M and the chart 
changings are defined on open subsets of the respective Frechet spaces. A C°°- 
atlas is an atlas with smooth chart changings. Two C°°-atlases are equivalent is 
their union is an C°°-atlas. A maximal C°°-atlas is called a C°°-structure on M 
(maximal is understood with respect to some carefully chosen universe of sets). A 
(smooth) manifold is a set together with a C°°-structure. A topological manifold 
is a manifold with a C'^-structure. If we are given a Banach space, then we can 
define C'^-manifolds by an C'^-structure. 

A smooth mapping f : M ~> N between smooth manifolds is defined in the 
canonical way, i.e. for any x G M there is a chart {V, v) with f{x) € V, a chart 
{U,u) of M with X G U and f{U) C V, such that v o f o u^^ is smooth. This is 
the case if and only if / o c is smooth for all smooth curves c : M ^ M , where the 
concept of a smooth curve is easily set upon. 

The final topology with respect to smooth curves or equivalently the final topol- 
ogy with respect to all inverses of chart mappings is the canonical topology of the 
smooth manifold. We assume manifolds to be smoothly HausdorfF (see the discus- 
sion in |T^, p. 265), i.e. the real valued smooth functions on M separate points. 
The product of smooth manifolds is defined canonically by building up the product 
of the atlases. 

A submanifold iV of a Frechet manifold M is given by a subset N C M, such that 
for each n £ N there is a chart (u, u{U)), a splitting E — E' x E" and u{U) — VxW 
with u{N) = V X {u{n)"}. By a splitting we shall always understand E' and E" 
as closed subspaces of E. 

A finite-dimensional manifold with boundary of dimension n is defined as or- 
dinary manifold except that we take open subsets in a halfspace {x S R" with 
a^n > 0}. For the notion (without surprises) of smooth mappings on such open sets 
see any textbook on differential geometry, for example The boundary {a; with 
Xn = 0} of the subspace models the boundary ON of the manifold TV, which is 
canonically a manifold without boundary of dimension n — 1 . A submanifold with 
boundary is given by the analogue submanifold-structure. 

We specialize to submanifolds with boundary of Frechet spaces - the case we 
are interested in. Let C -E be a submanifold with boundary, then the tangent 
space TnN at n £ N is given by the derivatives ■^\t=oc(t) of all smooth curves 
c :M. ^ N C E with c(0) = n. It is a vector space at any point n G N\dN and a half 
space at rt e dN. Nevertheless we denote by TnN the generated vector space and by 
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(T„A/')+ the canonically given halfspace. There is a natural structure of a manifold 
with boundary on the disjoint union UneNTnN , the tangent bundle. A vector field 
on iV is a smooth mapping on N associating to any point a tangent vector. A 
vector field X on C £^ is therefore simply a smooth mapping X : U ^ E, since 
the tangent space a.t x G U is E. The set of vector fields is denoted by X{N). The 
derivative Tf : TM TN of a smooth function / : M ^ TV is defined in the 
following way 

TJia,[c]) = {f{a),[foc]) 

where [c] denotes the equivalence class of curves c I'M. ^ M c E having the same 

derivative at the t = 0. A local diffeomorphism / : M ^ iV is an isomorphism on 
the respective tangent spaces. So the derivative mapping T^f : T^U Tf{x)V is 
given by T^f{v) = Df{x)v. 

Finite dimensional submanifolds (with boundary) of a Frcchct manifold M can 
locally be given by immersions. An immersion is a smooth mapping <j) : N ^ M 
with Ta(t> ■ TaN Tfj)(_a)M injective for all a G N and N a finite-dimensional 
manifold of dimension n. 

Lemma 4 (Submanifolds by Parametrization). Let N be a finite- dimensional ma- 
nifold (with boundary), E a Frechet space and (j) : N ^ E an immersion, then for 
any no £ N there is a small open neighborhood V such that (p{Vi) is a submanifold 
ofE. 

Proof. Given an immersion we shall construct a submanifold chart for the local 
image of N. We assume that N is an open subset of ffi.™ and - by translation 
- <;/)(no) = 0, since it is a local result. Given a linear basis ei,...,e„ of Tn^N 
we get linearly independent vectors Tna{4>){ei) =: fi £ E. We choose li,...,lm 
linearly independent linear functionals, such that li{fj) = 6ij and get a splitting 
E = E' X E" with dimE' = m via E" := kerZ^. The projection on the first 
variable pi induces a local diffeomorphism pi o (f> on a small open neighborhood 
V of no G TV by the classical inverse function theorem. The inverse is denoted by 
ijj : V C E' ^ V. Now we construct a new diffeomorphism 

r]{n, x") = {pi o (^(n), x" +P20 4>{n)) 

on y X W" , which is invertible by the above considerations: 

ri-\y',y") = my'U"-P2 0<t>my'm 

?7~^ defines a submanifold chart for ^(V^) since 

77-i(</.(n)) = (n,0) 

for n £ V hy definition. The proof for manifolds with boundary works in the 
same way except that the linear basis at a boundary point n € dN has to lie in 

{TnN) + . □ 

Two vector fields X G X{M), Y e 3C{N) are called /-related for a smooth 

map / : M — > TV if Tf o X = Y o f . We obtain a bounded linear mapping 
/* : ?i{N) — > X(M), the pull back, for a local diffeomorphism f : M ^ N hy the 
following formula 

{f*Y){x) = (T./)-i(y/(.)) 
for Y G X{N). The push forward for a diffeomorphism is defined by /* = (/~^)*- 
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Each vector field X admitting a local flow has a maximal local flow Fl produc- 
ing all integral curves on a connected set, furthermore X is FZj^-related to itself 
for any i, since 

T.Fl^{x) ■ X{x) = ±Flf{Flf{x))U^o = ^/lf+A^)U=o = X{Fl^{x)) 

Given a local semiflow, then the solutions of the differential equations are unique: 
Given an integral curve c we obtain 

^Fk-s{c{s)) = -X{Fk-s{c{s))) + X{Flt-s{c{s))) = 
as 

Fk_Ms))=c{t) 

for all s < t. 

Given two vector fields Y G X(M) and X e X(M) admitting a local flow, where 
M either denotes a finite-dimensional submanifold of a Frechet space or an open 
subset, we can define the Lie bracket, the most important notion of differential 
geometry: 

[X,Y] = ^^{FF,rYU^o 

Notice furthermore that for any local diffeomorphism (/) 

(t)*[X,Y] = [<t)*X,(t>*Y] 

so the pull back is a bounded Lie algebra homomorphism, since vector fields con- 
stitute a Lie algebra with the Lie bracket. For details in all possible directions see 

Finally we provide a formula for the Lie bracket on a Frechet space. Given two 
vector fields X and Y on an open subset U C E, where X admits a flow Fl^ , then 

[X,Y]{x)^j^{Fl^,rYix)U^o 

= jDFinFl^,{x)) ■Y{Fl^,{x))\,=, 

= DX{x) ■ Y{x) - D^Fl^{x){X{x),Y{x)) - DY{x) ■ X{x) 
= DX{x) ■ Y{x) - DY{x) ■ X{x) 

for X G U C E (see [|l^ and Due to this formula we can immediately conclude 

some properties of the submodule of Banach map vector fields: 

Lemma 5. Let U be an open set in a Frechet space E, then B{U) is a suhalgehra 
with respect to the Lie bracket. Let A be a bounded linear operator on E, then 
[A,B{U)] C B{U). Consequently the Lie algebra L{E) acts on B{U) by the Lie 
bracket. 

Proof. Given two Banach maps Pi andP2, DPi{x)-P2{x) = DQi{Ri{x)) ■ DRi{x) ■ 
P2{x) holds, which can be written as composition of DQi{v) ■ w for v,w G X and 
{Ri{x),DRi{x)-P2{x)) for a; G U. So the Lie bracket lies in B(f7). Given A G L{E), 
we see that APi {x) — DPi {x)-Ax is a Banach map by an obvious decomposition. □ 

Definition 4. Let E be a Frechet space, U an open subset. A distribution on U 
is a collection of vector subspaces D = {Dx}x£U of {TxU}x£U ■ A distribution is 
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called smooth if for any x £ U there is an open neighborhood x € V C U and n 
vector fields Xi,..., X„ with 

{X^{y),...,Xn{y)) = Dy 

for y E V. A distribution D on U is said to be involutive if for any two locally 
given vector fields X, Y with values in D the Lie bracket [X, Y] has values in D. A 
distribution is said to have constant rank ifdims.Dx is locally constant on U. 

We denote by (...) the generated vector space over the reals M. In the case of a 
smooth distribution D wc sometimes apply the notation D = {Xi, . . . ,X„) which 
means that is vector space generated by Xi{x),. . . ,Xn{x). We call this the 
linear span of Xi,. . . , X„. If the vector fields Xi,... , X„ are linearly independent 
in a open neighborhood of a point and if they span D, then we call these vector 
fields a local basis. 

Remark 7. Given a smooth distribution D on U such that the dimensions of 
are bounded by a fixed constant N . Let x G U be a point with maximal dimension 
Ux = dimR Dx, then there are n smooth vector fields X\,..., X„ and an open subset 
x GV CU such that 

{X,iy),...,X^{y))=Dy 

for y G V. At X there are linearly independent vectors, say Xi{x), Xn^{x). 
Choosing continuous linear functional e E' with li{Xj{x)) = 6ij, 

then the continuous m,apping I : U ^ L(W^^), x {li{Xj{x^|)) has range in the 
invertible matrices in a small neighborhood of x. Consequently in this neighborhood 
the dimension ofDy is at least n^- It follows by maximality ofux that it is constant. 
In particular a smooth distribution with constant rank admits locally a basis. 

The concept of weak foliations will be perfectly adapted to our purposes in HJM- 

theory: 

Definition 5. Given a distribution D on an open subset U C E, a tangential 
manifold (with boundary) for D is an immersion 4>pf : N ^ U such that C 
T„(f>N(T„N) for all n G N . We say that D admits a weak foliation of rank m if for 
any point xq € U there is an open subset V of {u G W" with Um > 0}, an open 
subset W" of a Frechet space E" and a smooth mapping ip : V x W" — > U such 
that 

i) There is a point in the parameter space such that xq = ^(uq, Xq); 

ii) T(^u.o,x")'>P '■ ^™ ^ E" E is an isomorphism; 

iii) %!}{., x") is an tangential manifold (with boundary) for D for any x" G W" ; 

iv) The tangent distribution given by all tangent spaces of the tangential manifolds 
has constant rank m. 

Remark 8. The immersed manifolds are called leafs. In the chart sometimes the 
notion plaque is applied for their local image. The dimension m of the tangen- 
tial manifolds will generically be greater than the dimension of of the given 
distribution D. 

Classically one is interested in the existence of weak foliations for a given dis- 
tribution of minimal dimension m. Therefore we shall need the following essential 
lemma: 
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Lemma 6. Let D be an involutive, smooth distribution of constant rank n on an 
open subset U of a Frechet space E. Let X and Y be vector fields with values in D 
and let X admit a local flow, then 

{Flfy{Y){x) G 

for X G U , where it is defined. 

Proof. Given a local basis Xi,...,Xn around xq, we have by involutivity that 
[X, X,] = fi^k- Remark that /f are smooth functions locally around xq: 

Given n linear independent functionals Im such that lm{Xj{xQ)) = 6m j, then 

n 

lU[X,Xi]{x)) = J2filmiXkix)) 

fc=l 

Since the matrix M{x) := {lm{Xk{x))) is invertible at Xq and has smooth entries, 
it is invertible on an open neighborhood of xq, and the inverse has smooth entries. 
The smooth inverse matrix applied to the left hand vector proves the smoothness 
of f^ . By involutivity we obtain 

|(FZf = -A(FZf_J*(X,)|.=o = 

= -^^iFl^X{Flfr{X,)U=o = 

= -[x.{Fifnx,)] = 
= -iFifr[x,x.] = 

n 

= -Y,f>Fl^{Fl^)\X,) 

k=l 

which is a linear equation with time-dependent coefficients gfit) := —f^{Flf{x)) 
for {Flf-)*{Xi){x) at any point x in that neighborhood of xq. Consequently pro- 
vided the initial values lie in the solution lies in D^, but {Fl^)* [Xi){x) = 
Xi{x). □ 

Now we can state and prove the Probenius theorem for weak foliations, remark 
that in this infinite-dimensional case the Frobenius condition is a sufficient condition 
for the existence of a weak foliation, since there are problems to solve ordinary 
differential equations. Nevertheless we shall see that there is a partial answer in 
the other direction. 

Theorem 8. Let D be an involutive, smooth distribution of constant rank n on an 
open subset U of a Frechet space E and assume that for any point Xq there is a 
basis of the distribution Xi, Xn of locally defined vector fields, where Xi, Xn-i 
admit local flows Flf* and Xn admits a local semiflow, then D admits a weak 
foliation of rank n. 

Proof. There are locally on some open set n linearly independent vector fields 
Xi, Xn generating each with local flows Fl-^^ for i = 1, n — 1 and a local 
semiflow Fl-^" . Without restriction we set xq = 0. We choose n linear functionals 
li, ...,ln such that li{Xj{0)) = Sij, so we get some splitting r : E ^ M.^ x E . We 
define the parametrization a{u,y) = Fl^^ o ... o Fl^j^{0,y) for (u, y) in an open 
subset around in R" x £ . We can calculate by looking at the partial derivatives 
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the tangent map applied to the canonical standard basis g|-|u '■— [s i— > u + sci] of 
(ruR") + , the 

Tua(^|u, g^lu)(^i. iFl^!)*-iFl^:-!)M{c^iu, y)) 

which lie by Lemma 6 in Da(u,y) and generate it in a small neighborhood by a 
dimension argument. It is essential that the first n — 1 vector fields admit a local 
flow. So we obtain a family of tangential manifolds for D. Each parametrization for 
fixed y defines locally a smooth submanifold with boundary a(ui, .., Un-i, 0, y). □ 



Remark 9. For details on Frohenius theorems in the classical setting see ]ll| ] . The 
phenomenon that there is no Frohenius chart is due to the fact that there is one vec- 
tor field among the vector fields Xi,...,Xn (generating the distribution D) admitting 
only a local semiflow. If all of them admitted flows, there would exist a Frohenius 
chart, which can he given hy a construction outlined in ]l9| . The non-existence 
of a Frohenius- chart means that the leafs cannot he parallelized, since they follow 
semiflows, which means that "gaps" hetween two leafs can appear. These "gaps" do 
not appear as interior point of a leaf. This is an infinite dimensional phenomenon, 
which does not appear in finite dimensions. 



Remark 10. Let D he a smooth distribution of constant rank n on an open subset 
U of a Frechet space E and assume that for any point xq there is a basis of the 
distribution X I, ...,Xn of locally defined vector fields, where Xi, Xn^i admit local 
flows Flf' and X^ admits a local semiflow. Assume furthermore that D admits a 
weak foliation of dimension n, then the distribution is involutive. This is easily seen 
hy applying (Fl^')* to Xj atxdU, which is by assumption an element of (since 
the flow restricts by uniqueness of integral curves to the tangential manifolds), so 
differentiation and smoothness yield the result by the formula of lemma 6. 



4. Existence of Finite dimensional Realizations 

Subsequently, we let assumptions (A1)-(A4) be in force. For simplicity we shall 
discuss only the time-homogeneous case (see Remark ^ for an outline of the time- 
inhomogeneous case). 

The vector fields {^,a^ , . . . ,a'^} induce two smooth distributions on D{A°°): 
their linear span D, and Dla, the linear span of all multiple Lie brackets of these 
vector fields. 



Recapturing the discussion in Section 1 we now can say that (1.1) has generically 
an FDR if D admits locally a weak foliation. We shall see below (Proposition ||) that 
the assumption that D admits a weak foliation on some open set V is equivalent to 
the assumption that Dj^j^ is an involutive, smooth distribution with constant rank 
on some open (in general smaller) subset U C V. 

Define the smooth map F := J^^i ■ ^ ^ D{A°°) by 

F^ (y) := 5(0^- (y)) - ^Dcj,^ (y) (^(</>^' (y))) . 

So that we can write fi{h) = Ah + T{£{h)). 

Let U denote an open set in D{A°°) in what follows. 
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Lemma 7. Let Xi, . . . ,Xk be linearly independent Banach-maps on U , for some 
A: G N. Then the set 

M = {h<^U\ fiih) G (Xiih), . . . , Xdih))} 
is nowhere dense in U . 

Proof. We argue by contradiction. Suppose there exists a set V C A/" which is open 
in D{A°°). For any h € V fl A/", there exist numbers ci{h), . . . , Cd{h) such that 

d 

Kh)^Y.''Ah)'^'{h). (4.1) 

i=i 

Since D{A*) is dense in H, we can find £,i, ■ ■ ■ ,£,d £ D{A*) such that the d x d- 
niatrix Mij{h) := ^i{a^ (h)) is smooth and invertible on V (otherwise we choose a 
smaller open subset V). Hence 







( um) \ 




= M-\h) 




V Cd{h) ) 







can be smoothly extended to V. Then (4.1) implies that A is a Banach-map on V. 
But this contradicts (A4), whence the claim. □ 

In what follows, we suppose that Dla has constant dimension ko 'EN on U . 

Lemma 8. We have ko > d and 

ii{h) i {(j\h),...,a'^{h)), VheU. (4.2) 

Moreover, for any ho G U there exists an open neighborhood V and Banach-maps 
X'^+'^,...,X^^-^ on V such that {/i, cr\ . . . , ct'', X'^+i, . . . , X'^^-i} is a basis for 
Dla on V . 

Proof. Assumption (A3) implies ku > d. By Lemma ^ there exists a nowhere 
dense set TV C ?7 such that /i(/i) ^ (o'^(/i), . . . , cj'^{h)), for aU /i e [7 \7V. Therefore 
ku > d. 

Now suppose ^{ho) G (o'^(ft-o), ■ ■ • , f^"^ (ho)) , for some ho G U. By the definition 
of Dla and Lemma there exist ku — d Banach-maps X'^+^j . . . ,X'^" on U such 
that 

DlaW = {a\h), a'ih), X^+\h), X'^- (/.)}, 

for h — ho, and hence for all h in a, neighborhood of ho, by continuity. But this 
implies that ii{h) lies in the span of Banach-maps, for all h d V. This contradicts 
Lemma 0, whence the claim follows. □ 

Remark 11. Lemma^ proves a conjecture in namely that any short rate real- 
ization is of dimension 2 (see Remark 7.1 therein). 

In view of Lemma ^ the minimal generic FDRs are of dimension d+1. The following 
result is an extension of ||, Proposition 7.4]. 

Proposition 1. Suppose ko = d + 1. Then there exist d linear independent (con- 
stant) vectors Ai, . . . , A^ G D{A'^) such that 

(tJ'(/i) G (Ai,...,Ad), VheU, j = l,...,d. 
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Hence necessarily we are in the (slightly modified) "deterministic direction vola- 
tility" case of ^ (see Section 6 therein). 

Proof. By (the proof of) Lemma |^ we have, for any k ^ 1, . . . , d, 

d 

[M,fT'^](/i)-^4XMf^-''(M, v/iec/. 

Fix hQ € U. As in the proof of Lemma we find an open neighborhood V of 
and £,1, ■ ■ ■ ,£,d S D{A*) such that the d x d-matrix Mijiy) :— £,i {<!>■' (y)) is smooth 
and invertible for all y e W := i{V). In view of (A2), the sets W and W :— 
{£, £ o A) {V) are open in B and B x B, respectively. A calculation shows that 

[ti,a']ih) = A\i{h),i{Ah)), 

where 

A'=(y, z) Aq^^iy) + OT(y)£(0'=(y)) - Dcj^'^iy) {z + l{r{y))) 

(linearity of £ is essential for obtaining this implicit dependence on h). Conse- 
quently, the functions 

d 

j^{y, z) M-^^y)^^ {^\y, z)) -.W -^^ 

i=l 

are smooth and satisfy 

d 

h\y,z)^Y.^'y{y,zW{y), 'i{y,z)eW'. (4.3) 



Differentiation of (4.3) with respect to z (which makes sense since W is open) 
yields 

d 

D<t>\y) = Y.D,j^{y,z)cj^{y), y{y,z) G W. 
i=i 

Arguing again by linear independence of {4>^, . . . , <j>'^}, we see that 

D,^^iy,z)^:(3^{y) 

are smooth functions of y only, with values in the dual space of B. Hence </> satisfies 
a linear differential equation on W, 

d 

We may assume that W is star-shaped with respect to yo = ({ho); that is, 

yn+t{y^yo)eW, ViG[0,l], VyeW. 

Otherwise we replace V by (^^{B^iyo)), where Be{yo) := {y G S | ||y - yoUs < e}, 
for some e > small enough. Let y G W and define ^p''{t) :— (f>''{yo + t{y — yo))- 
Then there exists an open interval /, which contains [0, 1] and such that 

d 

dti^'it) - J2 ^j^yo + ^iy - yo))iy - y^Wi^) 

V''(0) = /(yo), fc = l,...,d. 
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for t €z I. This differential equation has a unique solution, which is of the form 

d 

In particular, cj^^y) = 4(l)'/>'(2/o). 

We have thus shown that, for any y G W, there exists a d x d-matrix A'^{y) such 
that 

d 

</>'=(y) = 5]4(y)0^(yo). (4.4) 

By linear independence of the (p^s we conclude that {A^) is smooth and invertible 
on W . Now the proposition follows by a continuity argument. □ 

Let the assumptions of Proposition ^ be in force. We now shall see how an FDR 
looks. Note that we have 

Da\h)a^{h)G{Ki,...,Kd). e U. (4.5) 

Write 

v{h) Ah + anjAiih). 

By (^) it follows that Dla = ■ ■ ■ , ^d)LA = (^i ^ii • ■ • j ^d), where the latter 

equality can be deduced as in the proof of Lemma |8|. Also we see that 

[v,Kk]{h)^Dv{h)kk e (Ai,...,Ad), (4.6) 

for all h G U. Now let ho G U and A4 he a. leaf of the weak foliation of D in [/ 
through Hq (and hence a ((i+ l)-dimensional tangential manifold for (z^, Ai, . . . , A^)). 
As in the proof of Theorem we obtain a parametrization of A4 by 

d 

aiu,y)=Flu{ho) + ^ykAk, (u, y) e [0, e) x T/, (4.7) 

fc=i 

for some open set V cMf^ and e > 0, where Flu is the smooth semiflow induced by 
v. 

Though the FDRs are already nicely specified by the parametrization (4.7), the 
following geometric picture is worth noticing. Let A4 he a. leaf through Hq U as 
above. The band of d-dimensional affine submanifolds of A4, 

M{u) (Fluiho) + (Ai, . . . , Ad)) CiM, ue [0, e), 

forms a foliation of Ai. Using Taylor's formula we calculate, for h £ A4, 

v{h) = v{Flu{ho)) + f {Di^{Flu{ho) + s{h - FluihoMh - Fluiho))) ds 
Jo 

i^{Flu{hn)) + i){u,h). 

By ( ff.6[ ) we have i>{u, h) e (Ai, . . . , A^), for all h G N{u), and we already know 
th&t v{Flu{ho)) ^ (Ai, . . . , Ad). In other words, the vector-field restricted to 7V(w) 
splits into a component lying in (Ai, . . . , A^;) and a constant "outward pointing" 
component. 
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It remains to see whether more can be said about the functions Ai,...,Ad. 
We can write a'^ih) — X]^=i P'"' (^)^ji ^^r some smooth invertible matrix-function 
p{h) = Accordingly 

d 

i'{h)^Ah+ J2 a'^h)S{A„Aj), 
where a'^ (h) := {pp*Y^ {K) = ELi (^)- Now (U) imphes 

d d 

AKk+ (^a'^(MAfe)5(A„A,)=^7'^WA,. 

ij — 1 2—1 

Expressed as a point-wise equahty of functions this reads 

d d 

a,Afe(x) + T''''^{h)d^ (A,(a:)A,(x)) = ^7''(MA.(x), (4.8) 

ij = l 1=1 

where Ai(a;) := K{i)dv and T^^^^ [h) := {l/2)Da'^ {h)Ak. Integrating with 
respect to x yields 

d d 

d^Akix) = Ak{0) + ^7''(/i)A.(a:) - ^ r'=^*^"(/i)A,(x)A,(x), (4.9) 

i=l i.j = l 

and this has to hold for any h E U. Thus any h E U implies a system of Riccati 
equations for the functions Ai, . . . , Ad (which have to hold simultaneously). 
Suppose for the moment that the functions 

Ai, . . . , Ad and A^Aj, 1 < i < j < d, are hnearly independent. (4-10) 

Then ( ^ ) implies that "f^'^Qi) = 7*-'* and V^'"^^ [h) = F'^'*-' are constant on U. Choose 
continuous linear functionals vtj on D{A^) (they can be chosen from H) such that 

wy(Afe) = T^''^. 

Write 5'^{h) := a'^h) - v'^{h). Then we have 

D6'^h)Ak = 0. 

Hence d^^{h) = 6^mh]) is a function of the equivalence class [h] of h E U given by 
the equivalence relation f ^ g f — g & (Ai, . . . , A^). Thus 

a'3{h) = 6'^{[h])+v'^{h). (4.11) 



The parametrization (4.7) can always be globally extended for {u,y) e Ili>o x 
R'^ (though the image manifold can accumulate on itself). This does not mean, 
however, that a generic {d+ l)-dimensional realization for (1.1) exists globally. As 
an illustration we remark that a = (a*^ ) is by construction a strictly positive definite 
symmetric d x d-matrix. Thus (4.11) yields the constraints 

v''{h) + S"{[h])>0, i^l,...,d., 

which in general are not globally satisfied. For a concrete example see the following 
section. 
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Let US finally consider the particular case where p{h) = p (and hence , . . . ,a'^ 
are constant). Here ( ff.^ ) or ([f.S|), respectively, simplify to 

d 

d^Kk{x) ^Y.^''\h)K,{x). (4.12) 

i=l 

Since Ai, . . . , are linearly independent this implies that j^^{h) = 7*^' are con- 
stant. Equation (4.12) is a system of linear ODEs for the functions A,;. 

A necessary condition for constant dimension, kn — d+ 1, oi Dla on U is (4.2). 
This yields 



KM ^ (Ai,...,Ad), \/heU. 



(4.13) 



Integrating with respect to x it is easy to see that i^(h) G (Ai, . . . , A^) if and only 
if 



This again holds if and only ii h e O where 
O 



(4.14) 



{d d 
bl + Y.c,A,--J2a^^XA,\b, 
i=i *,j=i 



ci, 



Cd e 



(1 denotes the constant function l{x) = 1). We have dimDLA = d on O. Neverthe- 
less, O is a (d -t- l)-dimensional (affine) submanifold^ of £'(A°°). Hence, although 
Dla does not have full rank, d+ 1, globally on D{A°°) we obtain a global generic 
{d + l)-dimensional realization. That is, a global weak foliation of D{A°°) with 
leafs given by (4.7) and O, respectively. The global HJM-model is by construction 
Gaussian. 



Remark 12. The ODEs (^-^j md ((.ll), respectively, have already been derived 
and discussed in They are further analyzed in Q], where also the dynamics of 
the coordinate process is derived. 



Remark 13. We can interpret the coordinates of the invariant submanifolds eco- 
nomically. If we turn back to the construction of the charts as given in Theo- 
rem 8 we see that we have been choosing continuous linear functionals h, . . . ,ln- 
They can be chosen as linear combinations of point evaluations (we assumed them 
to be continuous) ev^^,... ,evx„ by Theorem 23.10 in [T2| , which tells that lin- 
ear combinations of point evaluations are dense m C°°(R>o,K) . For the proof 
look at the bounded multilinear mapping (eVx-^, . . . , eVx„) 1^ det^ev^. {Xj{y))) which 
cannot be identically zero, since then the extension of the multilinear mapping to 
(li, . . . ,ln) t— > Aei{li{Xj{y))) would be identically zero, too. However, this is im- 
possible due to linear independence of Xj{y). Consequently we can interpret the 
finite dimensional realization as Markovian process in "benchmarks" . Notice that 
even though "benchmark" are obtained by linear projections on the forward rate, 
the Markovian process of the finite dimensional realization is in general not the 
projection of the solution of the HJM- equation. 



^We assume here that A;, A^Aj G D(A°°). This requires a nice choice of 7*^' in (4.12). That 
is, the matrix (7'-') must have negative eigenvalues. 
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The necessity of the assumptions in Proposition ^ for the existence of generic 
FDRs is essentially given by 

Proposition 2. Suppose that D is a smooth distribution of constant rank d + 1 
on an open subset U . Assume furthermore that on U the distribution D admits 
a weak foliation of dimension m (by Definition |^ necessarily m > d + 1). Then 
diniK {Dla)x ^ ^'^^ there are points Xq G U such that D^a is a smooth involutive 
distribution of constant rank in a neighborhood of Xq . 

Proof. We repeat the argument of Remark ^ Take two vector fields of the canon- 
ical basis cr^, . . . , fj'*} of D . They admit a local (semi)flow, which restricts 
necessarily to the leafs of the weak foliation by assumption, so the formula of 
Lemma 6 tells that their Lie bracket lies in the tangent distribution of the weak 
foliation. By the special choice of the vector fields we know from section 2 that 
the Lie brackets admit local flows as Banach maps, so by uniqueness of the integral 
curve the flow restricts to the leafs. We can prove now by induction that Dla 
is an subdistribution of the tangent distribution of the weak foliation. Since the 
dimensions are globally bounded by m, it is a smooth distribution and therefore 
involutive around some points where the rank is maximal. □ 

The general case where D^a has locally constant rank ki) is sketched in the follow- 
ing proposition, which provides the full classification of FDRs in our framework. 
The detailed analysis and the geometric implications will be discussed elsewhere. 

Proposition 3. Suppose that Dj^a has locally constant dimension kjj and that the 
mappings 

{l,loA,...,loA'') : ^ S'^+i 

are open for all integers r > 0. Then there exist ko — ^ linearly independent vectors 
Ai, Afco-i e D{A°°) such that 

Dla = (m, Ai, Afc^_i) 

a^{h) e (Ai,...,Afe„_i) 

locally in h. 

Proof. We proceed as in the proof of Proposition 0. Fixing a point Hq, we can 
choose linearly independent Banach map vector fields Xi, Xk^^-i and an integer 
r > — 1 such that 

X'{h)^(t)'o{l,...,loA''){h) 

with (j)^ : V X <Z 5'"+^ D{A°°) smooth with V simply connected and 

Dla — {^^,Xi, ...,Xko~i) 

locally (a representation by constant vector fields is encoded with r = —1). We 
demand r to be the minimal integer with the above properties. That is, either 
r > and D^^cff ^ 0, for some i, or r = —1. We have to show that the latter is 
true. 

We argue by contradiction and suppose that r > 0. Notice that Xi, ....Xk^^i 
spans exactly the subdistribution of Dla generated by Banach map vector fields. 
Consequently we obtain [ii,X%h) = Y/jZi c'j{h)X^{h) locally in h, since the Lie 
bracket is a Banach map. We can write 

[f,,x'']{h) = A''m,-JoA^+\h)) 
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with 



-Dcl>'^{z^,...,zn 



/ /(r(^o)) \ 



due to the formula for the Lie bracket. Since we assumed that (/, / o A, I o ^4'"+^) 
is open we obtain 

A''{z°,...,z^+')= J2 lHz°,...,z^+')4r^{z°,...,zn 

i=i 

with smooth coefficients by the above arguments. Differentiation with respect 
to 2;''+^ and applying the differential to v G B yields 

kn-l 

D,.ct>\z\...,z^)-v= ^ {I3^{z\...,z^)-v)cl>'{z\...,z^), 
find a smooth curve c : IR — > B such that 



where D^r^^^z^, ^''+^) = Pj{z^, z'^). Given two points Zq and z'^ in B, we can 



c(0) = z5 and c(l) = z''. 

For fixed z'^ , z'^~^ we define V*^!*) •= '^''{z^^ ■■■■>z^~^,c{t)), a smooth curve into 
D{A°"), which satisfies the differential equation 

|V''=(i)= E fe(A...,^''-\c(i))|c(i))v^(t). 

j=i ^ 

This differential equation has a unique solution, namely there exist smooth curves 
A| such that 

Thus for t = 1 there are real numbers such that 



fen-l 



^'=(2°,...,^'-)= E s'<^'(^°,-,^'-\^5)- 

By smoothness and linear independence of the fields on both sides we conclude that 
there are smooth matrix-valued functions aj : V x B'^ ^ M., with existing smooth 
inverse, such that 

/(A-,^'')= E a^{z°,...,zn4>Hz°,...,z^-\z'o)- 

Hence <t)^{z^ , z'^'~^ , Zq), with frozen last variable Zq, are linearly independent and 
span on F X B"" the same subspace as the fields (j)''{z^, z^) do. Therefore we could 
have taken the with frozen last variable to span Dla- But this is a contradiction 
to the minimality of r. Hence r = — 1. That is, there exist constant vector fields 
Ai, Afc^_i around ho that span the sub distribution of Dla spanned by Banach 
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map vector fields. If we are given a basis of constant vector fields around any point 
/iQ, then we can extend the result to the whole domain of constant rank. This is 
the desired assertion. □ 

Remark 14. The time- dependent case is treated in a similar way: Generic FDRs 
for this invariance problem are essentially given by generic finite dimensional real- 
izations on the extended phase space D{A°°) = M x D{A°°). We look at the problem 
on the extended phase space with operators 

1 

fi{t,h) 


ffJ {t, h) 



m = 



Given an invariant manifold M in D[A°°), the projection pr2 '■ M x D{A°^) 
D{A°°) restricted to A4 is an immersion, since the second component of the vector 

fields on the extended phase space is spanned by 1. By invariance in D(A°°) we 
can conclude invariance in D{A°°). Given an invariant manifold A4 in D(A°°) we 
obtain a foliation by invariant manifolds inRx M for the homogenous problem on 
the extended phase space. 

On the extended phase space we can apply the Frobenius methods. See also |^ . 

Remark 15. We can replace D(A°°) in the above discussion by any Frechet space 
E that satisfies (H1)-(H3) (where H is to be replaced by E). The Banach-map 
(J : E ^ Eq has to be chosen according to (Al)-(A3). We assume that S{t) is a 
smooth semigroup on E and (A^) has to be replaced by 
(A4') A : E E is not a Banach map 
(see Lemma^. We may think of E ~ C°°(M>o;]R). Then all the above conclusions 
on the geometry of generic FDRs ( —weak foliations ) can be drawn. This illustrates 
that our analysis is essentially independent of the initial choice of H . The geometric 
problem even does not lead to more general solutions even on huge spaces of forward 
rates. 

5. The Svensson family as a leaf in a foliation 

A popular forward curve-fitting method is the Svensson family 

Gs{x, z) = zi + Z2e~'"^'-" + z^xe'^^"^ + 24X6^^'^'^. 

It is shown in j7| that the only non-trivial interest rate model that is consistent 
with the Svensson family is of the form 

rt = Zlg, + ■■■ + Zfgi, (5.1) 

where 

gi{x) = 1, g2{x) = e-""", gsix) = xe'""", 34(2;) = xe^^""^, 
for some fixed a > 0. Moreover, 

Zl^Zl^ Zf = Zle--\ Zt = Z^,e-'-* (4 > 0) 
and Z^ satisfies 

dZ^ = {Zf + Zf - aZf) dt + J^dWt. (5.2) 
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Here is a real-valued Brownian motion. 

We now shall find a generic local 2-dimensional realization that is of the form 
( 3.1) whenever rg = ^j9j with Z4 > 0. In view of (5^), a candidate for a is 

given, on U := {£ > 0}, by 



a{h) = y/ae{h)g2, 

where £ is some continuous linear functional on H (or even C(M>o, M)) with £{gi) — 
{^{92) = ^(53) — and ^(54) = 1 (notice that this is in full accordance with (4.11)). 
Straightforward calculations show, for h E U, 

H{h) = Ah + l{h)g2 - £{h)gl 



[^,cr](/i) = -ay/a£{h)g2 - —==g2. 
(the clue is that £ o = 0). Hence indeed dim{/i, u^la — 2 on \J . 
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